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the CMB and is uncovered in its bare form under the photon diffusion scale. Decoupling of the photons at last 
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oscillations, enhanced by the velocity overshoot effect, compete with CDM fluctuations in the present matter 
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of a growth-inhibiting radiation and baryon background. Combined with the acoustic contributions and 
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The sage sets aside obfuscation 

And is indifferent to baseness and honor 

The mass of men are all hustle-bustle 

The sage is slow and simple 
He combines myriad years 
Into a single purity 

Thus docs he treat the inyiiad things 
And thereby gathers them together 

-Chuang-tzu, 2 

I. Introduction 

Small scale fluctuations in the cosmic microwave background (CMB) and the matter density provide 
a unique opportunity to probe structure formation in the universe. As CMB anisotropy experiments reach 
toward smaller and smaller angles, the region of overlap with large scale structure measurements will increase 
dramatically. Since the CMB and matter power spectra encode information at very different epochs in the 
formation of structure, comparison of the two will provide a powerful consistency test for competing scenarios. 
Unfortunately, the simple relation between the two power spectrum at large scales (Sachs & Wolfe 1967) 
does not hold at small scales. To establish the general relation, we need to employ a more complete analysis 
of how gravitational instability affects the matter and CMB together. 

In previous papers (Hu & Sugiyama 1995a, b, hereafter HSa,b), wo developed a conceptually simple 
analytic description of CMB anisotropy formation. The central advance over previous analytic works {e.g. 
Doroshkevich, Zel'dovich & Sunyaev 1978) involved the influence of gravitational potential wells, established 
by the decoupled matter, on the acoustic oscillations in the photon-baryon fluid. In this paper, we refine 
the analysis for scales much smaller than the horizon at last scattering as is relevant for large scale structure 
calculations. We furthermore extend it to encompass the photon-baryon backreaction on the evolution of 
the cold dark matter and gravitational potential as well as baryonic decoupling and evolution. This makes 
it possible to extract the relation between CMB and matter fluctuations established in the early universe. 

Analytic solutions in terms of elementary functions can be constructed in the small scale limit and serve 
to illuminate the physical processes involved in a model-independent manner. Despite the simplicity of the 
final results, the study of small scale fiuctuations requires a rather technical exposition of perturbation theory. 
For this reason, we divide this paper into two components. The main text discusses results drawn from a 
series of appendices and illustrates the corresponding principles in the familiar context of adiabatic cold 
dark matter (CDM), adiabatic baryonic dark matter (BDM), and isocurvature BDM scenarios for structure 
formation. 

We begin in §2 with a discussion of the central approximations employed, i.e. the tight coupling limit 

for the photons and baryons and the external potential representation for metric fluctuations. Details of 
these methods and the justification of their use can be found in Appendix A. As established in Appendix 
B, all fluctuations are given a boost at horizon crossing due to the driving effects of gravitational infall 
and dilation. Since the gravitational potential subsequently decays, the driving effects disappear well after 
horizon crossing, leaving fluctuations to evolve in their natural source-free modes. For the photon-baryon 
system, these are acoustic oscillations whose zero point is displaced by baryon drag. As described in §3 
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and §4, these oscillations are frozen into the photon and baryon spectra at last scattering and the end of 
the Compton drag epoch respectively. Since these epochs are not equal, photon diffusion (Silk 1968) sets 
a different damping scale in the CMB and the baryons. This process is described by the acoustic visibility 
functions introduced in Appendix C. 

The baryonic oscillations may be hidden in CDM models by larger cold dark matter fluctuations. These 
are discussed in §5. The baryon and radiation background also suppresses CDM growth before the end of 
the Compton drag epoch. Appendix D treats this source-free evolution of the CDM component analytically. 
After the end of the Compton drag epoch, the CDM fluctuations provide potential wells into which the 
baryons fall. From the combined analysis of CMB, baryon and CDM fluctuations, we obtain accurate 
analytic expressions for matter and photon transfer functions in the small scale limit. This greatly improves 
upon the fitting functions for the matter power spectrum in the literature (Bardeen et al. 1986, hereafter 
BBKS, Peacock & Dodds 1994, Sugiyama 1995) in the case of a significant baryon fraction. Our form is 
constructed out of fits from Appendix E for parameters that depend on the detailed physics of recombination, 

1. e. the last scattering epoch, Compton drag epoch, photon damping length, and Silk damping length. 

2. Acoustic Approximation 

We previously developed an analjd;ic description of acoustic oscillations in the photons and baryons 
(HSa,HSb). As discussed in more detail in Appendix A, this approach is based on two simplifying assump- 
tions: that a pcrturbative expansion in the Compton scattering time between photons and electrons can be 
extended up to recombination and that the gravitational potential may be considered as an external field. 
Combined these two simplifications imply that photon pressure resists compression by infall into potential 
wells and sets up acoustic oscillations in the photon-baryon fluid which then damp by photon diffusion. 
However, both of these approximations would seem to be problematic for small-scale fluctuations. Here, the 
optical depth through a wavelength of the fluctuation at last scattering is small implying weak rather than 
tight coupling (see e.g. Kaiser 1984, Hu & White 1995). This issue is addressed more fully in Appendix 
C [below eqn. (C-3)] where it is shown that only acoustic contributions from the tight coupling regime are 
visible through the last scattering surface, i.e. those from a time slightly earlier than last scattering when 
the optical depth to Compton scattering was still high. 

On the other hand, the second issue poses a computational problem. Inside the horizon in the radiation 
dominated epoch, the acoustic oscillations in the photon-baryon density feed back into the evolution of the 
gravitational potential through the generalized Poisson equation. Hu & White (1996) take the approach 
of including the self-gravity effects consistently by separating them from the truly external potentials. For 
many cases, a simpler approach suffices. By considering the self-gravity of the fluid as contributing to an 
external potential, we retain the conceptual simplicity of HSa and HSb. Although a consistent evolution of 
the potential would seem to require a full numerical solution of the coupled equations, it is in practice not 
generally necessary. The key point is that the feedback effect on the potentials is not arbitrary. Since photon 
pressure prevents the growth of photon-baryon fluctuations and coUisionless damping eliminates neutrino 
contributions, the potential decreases to zero after horizon crossing in the radiation dominated epoch. We 
shall examine how this affects the acoustic oscillations more carefully in §3. 
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Fig. 1. Acoustic oscillations receive a boost at horizon crossing an due driving from 
gravitational potential decay. The perturbations then settle into a pure acoustic mode 
and are subsequently damped by photon diffusion. Together the potential driving and 
diffusion damping effects form the acoustic envelope. After diffusion damping has destroyed 
the acoustic oscillations, the underlying baryon drag effect becomes apparent. Since \E' is 
here constant for a » a^q, \R'^\ increases with time. This contribution itself is damped at 
last scattering a* by cancellation. Well after last scattering a ^ a* , isotropic perturbations 
coUisionlessly damp creating angular fluctuations in the CMB. The model here is adiabatic 
CDM, Oo = 1, /i = 0.5, Oft = 0.05. 

Due to Compton drag from the coupling, the baryons density fluctuations follow the photons as 
= |'^7 = 3G)o or 5b = 39o + S' and exhibit acoustic oscillations as well. Here overdots represent derivatives 
with respect to conformal time, Oq is the isotropic temperature perturbation in Newtonian gauge, and S 
is the constant photon-baryon entropy. On the other hand, the cold dark matter, if present, is decoupled 
from photons and suffers only the gravitational effects of the oscillating radiation. Again, since the potential 
created by the radiation merely oscillates and damps away after horizon crossing, the CDM density pertur- 
bations 5c experience a kick at horizon crossing only to settle into the pure logarithmically growing mode in 
the radiation dominated universe. 

3. Photons Fluctuations 

3.1 Driving Effects 

Before recombination, the photon-baryon system behaves as a damped, forced oscillator. As discussed 
in §2, the self-gravity of the acoustic oscillations contributes to the gravitational force on the oscillator at 
horizon crossing if the photons or baryons contribute significantly to the density of the universe at that epoch. 
Even though the detailed behavior at horizon crossing in the radiation dominated epoch is complicated by 
feedback effects, the end result is extremely simple: fluctuations experience a boost proportional to the 
potential at horizon crossing. As the first compression of the fluid is halted by photon pressure, the self- 
gravitational potential decays. This leaves the oscillator in a highly compressed state and enhances the 
amplitude of the subsequent oscillation. For adiabatic initial conditions, the amplitude of the potential at 
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horizon crossing is essentially the initial curvature fluctuation. For isocurvature initial conditions, it is simply 
related to the initial entropy perturbations. 

We quantify these statements in Appendix B and calculate the exact boost factor [see eqns. (B-5),(B- 
6)]. It represents an enhancement in temperature fluctuations of a factor of three over its initial value 
in the adiabatic case or more importantly a factor of 5/(1 + jgfu) over the gravitational redshift induced 
large-scale fluctuations in a scale-invariant model (Sachs & Wolfe 1967). Here the fraction of the radiation 
energy density contributed by the neutrinos is = Pv/{pv + P-y) — 0.405 for the standard thermal history 
with three massless species. Its presence accounts for neutrino anisotropic stress which provides the relation 
between the gravitational potential ^ and the curvature fluctuation $ [see eqn. (A-5)]. Note that neutrino 
temperature fluctuations are boosted by a similar factor (Hu et al. 1995). 

Three features are worth emphasizing: 

(1) Potential decay enhances temperature perturbations. 

(2) Enhancement only occurs for scales that cross the horizon before equality. 

(3) It occurs at horizon crossing causing phase-coherent oscillations in the wavenumber k. 

Together the first two points imply that as a function of fc, the acoustic amplitude will increase through k^q, 
the horizon crossing wavenumber at equality, forming a potential envelope. By delaying equality, one moves 
the transitional regime to larger scales and enhances the lower k oscillations. This explains the increasing 
prominence of the first few acoustic peaks as the matter content Qoh^ is lowered or the number of relativistic 
species raised in adiabatic models. It is interesting to note that the effect of equality on the CMB and matter 
power spectrum are anticorrelated providing a powerful consistency test for its redshift Zgq. As we shall see 
in §3.2 though, a complication arises since high-fc oscillations are damped by photon diffusion. 

The third point follows because the gravitational driving term is effective at sound horizon crossing 
r] ~ 1/fccs, thus mimicing a driving frequency of u) ^ kCg- Here rj — J dt/a is the conformal time and we 
take c = 1 throughout. Since the natural frequency of the oscillation is related to the soimd speed c,, as 
w = kcg, the two scale in the same manner. Thus the horizon crossing effect is timed to produce a coherent 
oscillation in k. Specifically, adiabatic and isocurvature initial conditions yield cosine and sine temperature 
oscillations respectively. 

3.2. Damping Effects 

The photons and baryons are not in fact perfectly coupled leading to diffusion damping. The coupling 
strength is quantified by the Compton optical depth r. Correspondingly, the mean free path of the photons 
in the baryons is given by f~^. Explicitly, f = XeUeCTTa where is the electron ionization fraction, rig 
is the electron number density, and or is the Thomson cross section. As the photons random walk across 
a wavelength of the perturbation, temperature fluctuations collisionally damp. More specifically, diffusion 
generates viscosity or anisotropic stress in the photon-baryon fluid and causes heat conduction across the 
wavelength (Weinberg 1972). Both of these processes damp fluctuations. To order of magnitude, the damping 
length is the random walk distance yjrjjf. It exceeds the wavelength of a fluctuation when the optical depth 
through a single wavelength f/k = krj. If the diffusion scale is well under the horizon, krj ^ 1 so that f/k 
is still high at crossover. The photons and baryons are thus still strongly coupled and the damping may be 
calculated under the tight coupling approximation. 
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A quantitative treatment of diffusion damping is given in Appendix A. Subtle effects such as the angular 

and polarization dependence of Compton scattering slightly enhance the generation of viscosity and thus 
damping in the radiation dominated universe (Kaiser 1983). The end result of the calculation is a wavemimber 
koiv) by which acoustic fluctuations arc damped as cxp[— (fc/fc/j)^] [sec eqn. (A-14)]. Combined with the 
potential envelope from the horizon crossing boost of §3.1, this completes the acoustic envelope shown in 
Fig. 1. Remaining after diffusion damping is the acoustic offset of 8o + ^I/ = —R'i, where R = 2>pb/Ap^ = 
31.5f26/i^9^7(2:/10^)"^ Here 62.7 = To/2.7K is the scaled present temperature of the CMB, Ob is the 
fraction of critical density contributed by the baryons, and the Hubble constant is Hq = lOO/i km s~^ 
Mpc~^. We add the Newtonian potential ^ to the temperature perturbation ©o to remove the effect of 
gravitational redshift on the photons (see Appendix A). The term represents the baryon drag effect on 
the photons and is analogous to the Silk mechanism (Silk 1968) with the roles of the photons and baryons 
interchanged. Infalling baryons drag the photons into potential wells leading to a displacement of the zero 
point of the acoustic oscillation (HSa) . Thus it is responsible for the alternating peak heights and amplitude 
enhancement of intermediate scale acoustic oscillations (see Appendix A, Fig. 8a). The zero-point shift 
remains even after diffusion damping has eliminated the oscillations themselves. 

3.3. Decoupling 

As the total optical depth to the present drops below unity r(2;*) = 1, last scattering of the CMB photons 
freezes the acoustic oscillations into the spectrum. The optical depth drops rapidly as neutral hydrogen forms 
so that last scattering and recombination approximately coincide in the absence of subsequent reionization. 
This epoch is nearly independent of cosmological parameters such as the matter and baryon content Cloh^ 
and fift/i^ (Peebles 1968, Jones & Wyse 1985). However variations at the 10% level do occur across the full 
range of parameters. In Appendix E [eqn. (E-1)], we present an extremely accurate analytic fit to the last 
scattering epoch. 

The phase of the acoustic oscillation is frozen in at the value fcrs(?7,), where (HSa) 



is the sound horizon with R^g = R{aeq) and Ugq = 2.35 x 10 ^{Qoh'^) ^(1 — /y) ^©It- Its variation with k 
produces an oscillatory pattern in the CMB temperature with extrema at scales (HSb) 



The physical scale of the peaks k-yj is projected onto an angular scale on the sky and provides a sensitive 
angular diameter distance test for curvature in the universe. 

To determine the amplitude of the fluctuations, this instantaneous decoupling approximation must 
be modified to account for diffusion damping through recombination. The differential visibility function 
V-y = fe""^ expresses the probability that a photon last scattered within drj of 77. This function describes 
how fluctuations at recombination are frozen into the CMB. Equally useful for our purposes is the acoustic 
visibility function = exp[— (fc/A;£))^]. By including the intrinsic damping behavior of the fluctuations, 
it describes how acoustic oscillations are frozen into the CMB (see Fig. 2). Due to the growth of the diffusion 
length through recombination, this function is weighted to slightly earlier times than V-y. At small scales. 
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Fig. 2. Visibility functions for the photons and baryons. In the undamped fc — > limit, the 

photon acoustic visibility and the Compton visibility are equivalent = and the baryon 
acoustic visibility equals the drag visibility, Vb = V^. If ^bh? < 0.03, Vb peaks at later times 
than V-y, i.e. r/d > r]t,. For small scales, the acoustic visibilities, which weight the acoustic 
contributions from the tight coupling regime, arc increasingly suppressed at later times by 
damping and the width and amplitude of the acoustic visibilities decrease as k increases. 

only the small fraction of photons last which last scattered before z, , and hence in the tight coupling regime, 
retain acoustic fluctuations. This leads to a sharp decrease in acoustic fluctuations with k. More specifically, 
the damping envelope is given by (HSa) 



V^{k)= / dri%{'n,k)= / dr]V^{ri)e 
Jo Jo 



-[k/kniv)f 



(3) 



i.e. a near exponential damping in k. Through the first decade of the drop, it is well approximated by the 
form 

P^(A:) ~ e-['=/'=°^l"'\ (4) 

The damping angle is given by the same projection conversion as the acoustic peaks £d = kD'yr${T}^). 
Accurate fitting functions for ko-y and are given in Appendix E, cqns. (E-4) (E-7). Note that ku^ is 
the wavenumbcr at which diffusion suppresses the fluctuation by e^^. It is interesting to note that k^l^ is 
intimately related to the width of the Compton visibility function. This is because the thickness of the last 
scattering surface is by definition the diffusion length at last scattering. 

Combining the damping envelope with the intrinsic amplitude of acoustic oscillations discussed in §2 
and Appendix B, we obtain the transfer function at small scales. This function represents the growth to 
the present from the initial curvature or entropy perturbation. In Fig. 3, we plot examples for the adiabatic 
BDM, isocurvature EDM, and adiabatic CDM models with standard recombination. Note that in more 
typical isocurvature BDM models (Peebles 1987a,b), reionization may occur soon after recombination. In 
this limit, the damping length continues to grow until it reaches the horizon at the new last scattering surface 
and destroys all oscillations in the photons. 

The corresponding anisotropy can be obtained by choosing an initial curvature and entropy spectrum 
$(0, k) and 5(0, k) and employing the free streaming solution (see HSa eq. 12) for the monopole and dipole 
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Fig. 3. Photon transfer function. Plotted is the present rms photon temperature fluctuation 

relative to the initial curvature fluctuation $(0, /c) for adiabatic fluctuations and the initial 
entropy fluctuation S(0, k) for isocurvature fluctuations assuming standard recombination. 
The intrinsic acoustic amplitude is approximated by the analysis in Appendix B for scales 
well inside the horizon at equality k 3> k^q. The damping is well approximated by the tight 
coupling expansion for high Q.})h^ and is slightly overestimated for low values. Below the 
damping tail, the baryon drag offset clearly appears in (b) and (c) where the gravitational 
potential is not dominated by acoustic density fluctuations. An analytic treatment of this 
effect is given in Appendix A and C. 

contributions to the rms given in Appendix B, eqn. (B-7). Note that for low Q.bh? models, the damping length 
is somewhat overestimated by the tight coupling approximation. This is not surprising since as Vli,h^ — > 0, 
the photon mean free path approaches the horizon. In this case, the diffusion length passes the wavelength of 
the fluctuation when the optical depth through a wavelength is near unity and the tight coupling expansion 
of Appendix A breaks down. The damping length is oiierestimatcd because the photons essentially free 
stream across the wavelength and do not suffer coUisional damping. A phenomenological correction for this 
effect is given in Appendix E. 

In summary, the general features of the CMB fluctuation spectrum are 

(1) Acoustic peaks at harmonics of the angle the sound horizon subtends at last scattering, sensitive to 
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Fig. 4. Baryon fluctuation time evolution. The baryon density fluctuation Sb foUows the 
photons before the drag epoch yielding a simple oscillatory form for an <C a <C a^. The 
Silk damping length is given by the diffusion length at the drag epoch. The portion of the 
baryon fluctuations that enter the growing mode is dominated by the velocity perturbation at 
the drag epoch due to the velocity overshoot effect (see also Fig. 11). Since a flat A model 
is chosen here, the growth rate is slowed by the rapid expansion for a> a\ = [flo/fl^y/^. 

(a) The curvature, 

(b) The adiabatic or isocurvature nature of the initial fluctuations; 

(2) An acoustic envelope for the oscillation amplitude sensitive to 

(a) The matter-radiation ratio, e.g. the adiabatic acoustic envelope varies from to as one passes 
through the equality scale, 

(b) The baryon-photon ratio due to the a modulation of from baryon drag, 

(c) The thermal history from the exponential diffusion- damping tail, 
though the exact nature will depend on details of the model. 

4. Baryon Fluctuations 

In the tight coupling limit, the baryons participate in the acoustic oscillations of the photons. Near 
recombination, they decouple from the photons though not exactly at last scattering. Scattering represents 
an exchange of momentum between the two fluids and seeks to equalize their bulk velocities Vb and Vj. 
However the two momentum densities (p-y + Pj)Vj = ^p-yV^ and {pb + Pb)Vb pbVb are not equal. Thus 
momentum conservation requires that the rate of change of the baryon velocity due to Compton drag is scaled 
by a factor of = ^p^/pb compared with the photon case: i.e. oc fa = t/R. The explicit expression for 
the baryon momentum conservation or Euler equation is given in Appendix A and is used in Appendix C to 
make the qualitative statements here rigorous. The form of the coupling suggests that we can define a drag 
depth Tdirjd) = J^^^ Tddrj. Below drag depth Td{rjd) = 1, the baryons dynamically decouple from the photons. 
For the standard recombination scenario, this occurs near recombination but at a different value than last 
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Fig. 5. Matter transfer function. The analytic estimates of the intrinsic acoustic amphtude 

is a good approximation for k ^ kgq. The Silk damping scale is adequately approximated 
although its value is underestimated by ~ 10%. For isocurvature BDM and adiabatic CDM, 
the acoustic contributions do not dominate the small scale fluctuations. We have added in 
the contributions from the initial entropy fluctuations and the cold dark matter potentials 
according to the analytic treatment of Appendix E. 

scattering. Analytic fitting formulae for zj, are given in Appendix E, equation (E-2). Since recombination 
occurs around z = 10^ for all models, the end of the drag epoch precedes last scattering if Qbh"^ ^ 0.03. 

The acoustic fluctuations in the baryons are frozen in at the drag epoch rather than at last scattering. 
Furthermore, unlike for the CMB, it is not the acoustic density fluctuation that forms peaks in observable 
spoctriim today but rather the acoustic velocity. This is because the baryon fluctuations continue to evolve. 
In Appendix C, we give the exact matching conditions at Zd onto the growing and decaying modes of pres- 
sureless perturbation theory. This yields an accurate description for the subsequent evolution of the baryonic 
fluctuations in the presence of a background radiation energy density and cold dark matter. Qualitatively, 
the acoustic velocity at the drag epoch dominates over the acoustic density for the growing mode of fluctu- 
ations due to the velocity overshoot effect (Sunyaev & Zel'dovich 1970, Press & Vishniac 1980). The former 
moves matter and produces clumping in the baryon density. Since expansion damps peculiar velocities, this 
lasts for approximately an expansion time Thus only scales smaller than the horizon krjd ^ 1 experience 
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a boost due to the velocity at release. We show an example in Fig. 4, where the fc-mode is chosen to be near a 
zero point of the acoustic density oscillation at 77^;. The rapid regeneration of density fluctuations via velocity 
overshoot is due to the fact that zeros of the density oscillation are maxima of the velocity oscillation. The 
peaks in the matter power spectrum due to baryonic acoustic oscillations therefore occur at 

ku r (n.) - / ~ ^Z^)'^' adiabatic . . 

kbjrsind) - I .^^ isocurvature 

and are roughly n/2 out of phase with the corresponding CMB fluctuations. 

To obtain the amplitude of the acoustic fluctuations, we must also consider damping effects. Photon 
diffusion in the tight coupling regime damps baryon fluctuations as well due to Compton drag, i.e. via the 
Silk mechanism (Silk 1968). Analogous to the photon case, we can construct the drag visibility function Vt 
out of the drag optical depth t^. The acoustic visibility function then becomes Vi, = Vb exp[— (fc/fci))^] (see 
Fig. 2). Similarly, the net damping as a function of scale is described by 

Vb{k)= / dr]Vb{v,k)= / d7?V6(r?)e-t'=/'=°('''l ~ e-(*=/'=«)'"^ (6) 
Jo Jo 

where the approximation is valid through the first decade of damping. Analytic fitting formula for ks and 
ms are given in Appendix E eqns. (E-9) (E-10). As is the case with the photons, the latter accounts for the 
width of the visibility function and is almost independent of cosmological parameters. Note that the Silk 
damping length is not the same as the photon damping length despite the underlying similarity in cause. 
The small difference between last scattering and the drag epoch can alter it significantly due to the rapid 
change in ko around recombination (see Fig. 2). Moreover, the two scale differently with the baryon and 
matter content. 

Together with the horizon crossing boost from Appendix B, this defines the contribution to the matter 
transfer function of the acoustic oscillations (sec Fig. 5). Explicit expressions for the three scenarios shown 
arc given in Appendix D eqns. (D-22) (D-23). In the isocurvature EDM models, the transfer function at small 
scales is dominated by the initial entropy fluctuation 5(0, A;). Furthermore, unlike the photon oscillations, 
baryon oscillations may survive early reionization if it occurs more than an expansion time after the drag 
epoch. In this case, the baryonic oscillations are subsequently surrounded by a homogeneous and isotropic 
distribution of photons. They then represent entropy perturbations and are not damped by further photon 
diffusion. Even in the rather unphysical event of near instantaneous reionization, baryonic oscillations may 
survive if reionization is accompanied by the formation of a significant fraction of compact baryonic objects 
(see e.g. Gncdin & Ostrikcr 1990). 

If the model contains cold dark matter, baryons suffer an additional effect. After the drag epoch, 
they fall into potential wells established by the CDM. If the CDM to baryon ratio is high, this effect will 
dominate over the velocity overshoot of the acoustic oscillations. To quantify this effect, we need to consider 
the evolution of CDM fluctuations. 
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Fig. 6. CDM and matter fluctuation time evolution. The cold dark matter fluctuations 

are constant outside the horizon scale and experience a boost into a logarithmically growing 
mode at horizon crossing in the radiation dominated epoch. Between equality and the drag 
epoch, the presence of baryons suppress the growth rate of CDM fluctuations [see cqn. (7)]. 
By lowering fio^^j this region of suppression can be reduced for fixed fi^/rio, c/. (a) and 
(b). After the drag epoch, we plot — [SbPb + ^cpc)/{pb + Pc) rather than dc since the 
two components thereafter contribute similarly to the total growth. Baryons also lower the 
contribution of Sc to 6m at the drag epoch. 

5. Cold Dark Matter Fluctuations 

Let us begin with the evolution of CDM fluctuations before the drag epoch. As shown in §2 and 
Appendix B2, CDM fluctuations are given a kick at horizon crossing that sends them into a logarithmicly 
growing mode. As the universe becomes matter dominated, this stimulates power-law growing and decaying 
modes 5c oc a^. If CDM dominates the non-relativistic matter, p — {1, —3/2}. However if the baryon fraction 
is significant, the power law is modified. To first order in Qb/^o, 



, 3 rife 3 
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CDM growth is thus inhibited by the presence of baryons. We give the exact solution to the evolution 
equation from horizon crossing to Zd in Appendix D. Because of the complexity of these expressions, it 
is also useful and instructive to obtain approximate scaling relations. If Ueq and Qb/^o ^ 1, 

the main effect is an amplitude reduction of the CDM density perturbation 6c{rid,k) by approximately 
{ad/aeq)~°-^^'-^^° ^ {24floh'^)-°-^^f-/^° (see also Fig. 10). As floh'^ is lowered, the drag epoch recedes into 
the radiation domination epoch and the regime where the growth rate is affected Ueq < a < ad vanishes. 

At the drag epoch, the baryons are released from the photons and behave dynamically as if they were 
CDM. Since baryonic infall into CDM wells subsequently contributes to the self-gravity of the matter, the 
growth rates again become p = {1, —3/2} regardless of the baryon fraction. We thus follow the total non- 
relativistic matter perturbation dm after the drag epoch. The relative contribution of the CDM fluctuations 
at the drag epoch to the total non-relativistic matter fluctuations Sm scales as 1 — i^b/^o- A good flt to the 
net suppression is given by the form 



a 1 



fife 



X lim Sr, 



(8) 
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where a ~ {AlQ.oh?)-^-^'^^''/^'' for low Q.b/O.o < 0.5 and high Q.Qh? » 0.03. The change in the coefficients 
from the naive scaling relation is due to detailed matching of growing and decaying modes [see Appendix D 
eqn. (D-18)]. A highly accurate fitting formula for a in the general case is given in Appendix E, cqn. (E-12). 
For the ^ limit, an exact expression in terms of elementary functions is given in Appendix D which 
improves the 10% accuracy of the standard BBKS fitting formula to better than the 1% level at small scales. 

In Fig. 6, we show the resulting evolution of a scale under the Silk damping length for which the baryon 
fluctuations at the drag epoch are negligible. To extend the scaling of equation (8) to larger scales for Fig. 5, 
we have employed a generalization of the BBKS fitting function given in Appendix D. Notice that 

(1) the change in the growth rate between equality and the drag epoch, and 

(2) the fractional contribution of 5c to 5m at the drag epoch, 

both play a significant role in suppressing the final amplitude of matter fluctuations. Combined with the 
acoustic contributions from §4, this completes the matter transfer function in CDM models. 

We can now address the question of when acoustic oscillations are prominent in the matter transfer 
function. The main effect is simply due to the density ratio Pb/ pc = ^b/^c- However, the acoustic and 
cold dark matter contributions have a different dependence on scale. Relative to the cold dark matter, 
acoustic contributions scale as {krid)'Db due to the velocity overshoot and diffusion damping factors. Since 
Vb encorporates an exponential cut off at the Silk scale ks and velocity overshoot weights the spectrum 
toward small scales, acoustic contributions will be most visible just above the Silk scale. The relative 
contribution to the matter transfer function will therefore scale as ksr}d- Acoustic contributions increase in 
prominence if the Silk scale is small compared with the horizon at the drag epoch, i.e. in the high Q.bh? case. 
By including the suppression factor from equation (8) and numerical factors from Appendix B and D, the 
maximum ratio of the acoustic amplitude to the CDM contribution in the transfer function scales crudely 
as [see Appendix D, eqn. (D-31)] 

0.4fts^(f^o/i')"'(l + 24rJo/i' + 32r!f,/i2)-3/4a-i, (9) 

* 'c 

where ks, here in Mpc~^, and a are given explicitly in Appendix E. This relation encorporates: 

(1) Baryon acoustic oscillations; 

(2) Baryon decoupling; 

(3) Silk damping; 

(4) Velocity overshoot; 

(5) Baryon gravitational instability; 

(6) Baryon infall into CDM wells whose depth depends on 

(a) Acoustic feedback into the potential, 

(b) CDM self-gravity vs. the expansion rate, 
and combines them in a consistent manner. 
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6. Discussion 

We have established a framework for treating small scale fluctuations in the realistic case of a coupled 
multifluid system. In the limit that fluctuations crossed the horizon in the radiation dominated epoch, 
closed-form analytic solutions arc available. The fundamental elements uncovered by this approach arc the 
boost at horizon crossing due to infall and dilation effects from potential decay, the source-free solution of 
the component evolution equations, the baryon drag on the photons, and the Compton drag on the baryons. 
Together they establish general consistency relations between the matter and the radiation power spectra as 
well as expose their sensitivity to changes in the background model. 

Unlike the case of the matter fluctuations, decay in the potential results in an amplification of acoustic 
oscillations in the CMB. These opposing manifestations of the same physical effect provide both a measure 
and a powerful consistency check on the redshift of equality. Baryon drag displaces the zero point of the 
acoustic oscillations which remains as a temperature shift even after the oscillations have damped by photon 
diffusion. Last scattering marks the end of the baryon drag epoch at which the acoustic oscillations with 
their characteristic diffusion damping scale are frozen into the photon spectrum. Correspondingly, at the 
end of the Compton drag epoch, the Silk-damped baryonic acoustic oscillations are frozen into the matter 
spectrum. We have provided convenient analytic fitting formulae for these quantities as a function of the 
matter and baryon content. These may be useful for the extraction of cosmological information from the 
CMB and matter power spectrum once observations become available. 

The photon-baryon system also affects the growth of CDM fluctuations. It first stimulates growth 
through the horizon crossing boost. Subsequently, it affects the balance of the growth inhibiting expansion 
to the self-gravity of the CDM. We have obtained an analytic solution for adiabatic initial conditions and the 
exact general solution for the two effects respectively. To simplify these expressions, we have also provided 
an accurate fitting form for the transfer function in terms of elementary functions. Growth suppression due 
to the presence of baryons has implications for the first generation of structure. The expressions derived 
here remain valid for linear perturbations to an arbitrarily small scale where direct numerical calculations 
are impractical. 

Baryonic acoustic oscillations are of course not prominent in presently popular models where 0(,/f2o ^ 1 
due to the extra growth of the CDM fluctuations between horizon crossing and the drag epoch. However, 
they serve as a useful complement to their CMB counterpart if either floh^ < 0.05 or big bang nucleosynthesis 
constraints arc too stringent (Gnedin & Ostrikcr 1990). Indeed there arc tentative indications from cluster 
inventories that the baryon fraction may be as high as 15% (White et al. 1993). The presence or absence 
of acoustic oscillations in the observations of the CMB and large scale structure will in the future provide a 
robust distinction between general classes of scenarios: 

(1) Oscillations in CMB and matter power spectra: standard recombination with O;, > Ocl 

(2) Oscillations in CMB alone: standard recombination with fib <^ Qc', 

(3) Oscillations in matter power spectra alone: early reionization with f2(, > Oc- 

Large scale structure observations already suggest there are no dramatic oscillations in the matter power 
spectrum as would be the case for (1) and (3) if fib ^ flc (Peacock & Dodds 1994). However low amplitude 
oscillations, as might be expected if fib ~ ^c, remain possible. Of course, the exact form that these 
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oscillations will take in the observations depends on issues such as redshift space distortions and non-linear 
corrections. 

If oscillations are discovered in neither spectra, the most natural conclusion is our universe has fij ^ VLc 
and suffered early rcionization. However other possibilities include the formation of perturbations after 
recombination, rcionization within an expansion time after the drag epoch, equal or random stimulation of 
adiabatic and isocurvature mode acoustic fluctuations.^ These scenarios can be distinguished by measuring 
the location of the damping scale. All scenarios obeying standard recombination, regardless of the presence 
of actual peak-like structures, follows the scalings for the damping length discussed here. In all reionized 
models, the horizon at last scattering and at the drag epoch marks the damping scale for the photons and 
baryons respectively. 

If acoustic oscillations are discovered in both the CMB and large scale structure power spectra, we 
will possess a strong consistency test for the dynamics of the expansion, i.e. a combination of the matter 
content, curvature, and cosmological constant, as well as the adiabatic or isocurvature nature of the initial 
fluctuations. Furthermore, the two contain complementary information on several speciflc fundamental 
cosmological parameters. The scale of the peaks in the matter power spectrum are mainly determined by 
the matter content Oq^^ whereas the angular scale of the CMB peaks is mostly sensitive to the curvature 
1 — flo — ^^A- The two damping lengths also probe diflferent combinations of Q.Qh? and Ofc/i^. The ratio of the 
peak heights to the underlying CDM contribution in the matter power spectrum probes Ob/Qo- Furthermore, 
by comparing similar scales, dependence on the initial power spectrum can be eliminated providing a clean 
test of the whole gravitational instability paradigm. 
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Appendix A: Tight Coupling Limit 

A.l. Evolution Equations 

The Fourier transform of the Newtonian temperature fluctuation can be broken up into Lcgendre mo- 
ments, related to the direction cosines of the photon momenta 7^, 6(77, k, 7) = ^^{—iYOgPiik ■ 7). The 
evolution equation for these moments is given by the Boltzmann hierarchy (Bond & Efstathiou 1984, HSb) 



eo = -^ei-$. 



Gi = A: 
62 = A; 
ee = k 



60 + * --82 
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(A-1) 



10 ^ 2 
- fee, {£ > 2) 



if the ratio of the wavelength to the curvature scale is much smaller than the angle considered, i.e. k ^ t^/—K 
where K = — iJg (1 — fio — Oa)- Recall that overdots are derivatives with respect to conformal time 77 = / dt/a. 
Here 8^ is the CMB temperature perturbation in the Stokes parameter Q. It accounts for polarization 
generated by Compton scattering of anisotropic radiation. The metric fluctuations in the mode are given 
by 500 = ~(1 + 2^'y) and gij = a?{l + 2$y)7ij, where 7ij is the three metric on a surface of constant 
curvature and y is a plane wave e'"*'^ in flat space or more generally a fc-eigenfunction of the Laplacian. The 
presence of the curvature perturbation $ in the monopole equation represents the dilation effect. The form 
of the metric shows that it has the same origin as the photon redshift with the expansion. The gravitational 
potential 5* in the dipole or velocity equation accounts for gravitational infall or redshift. 

The tight coupling approximation assumes that the Compton scattering rate f is sufficiently rapid to 
equilibrate changes in the photon-baryon fluid. It is an expansion in the Compton scattering time f^^, or 
more specifically the inverse of the optical depth through a wavelength f/k and through a period of the 
oscillation f/oj = f/kcg > f/k, where 

Cs = , ^ (A-2) 
^3(1 + ^) 

is the photon-baryon sound speed with K = jp^/p^. To first order, only the £ = monopole (with density 
fiuctuation 6^ = 48o) and £ = 1 dipole (with bulk velocity = 61) survive and one obtains the forced 
oscillator equation for acoustic waves in the photon-baryon fluid (HSa).t To second order, the acoustic 
oscillations of the monopole and dipole are damped due to the imperfect coupling between the photons 
and baryons. Photon diffusion creates heat conduction through 81 — Vb and shear viscosity through 82 
(Weinberg 1972, Bond 1995). 



t In HSa and HSb, we employed a hybrid gauge or "gauge invariant" representation of density fluctuations 
for computational convenience. Since there are no benefits of this choice below the horizon, we work entirely 
in the Newtonian gauge in this paper. Only the definition of density fluctuations is affected: the total matter 
gauge Ax = 5x + 3^(1 +Px/px)VT/k where X represents any of the particle components. 
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To close these equations, we need the continuity and Euler equations for the baryons 

5b = -kVb - 3$, 
Vb 



--Vb + k^ 
a 



f(ei - Vb)/R, 

the polarization hierarchy equations for the CMB (Bond & Efstathiou 1984, Kosowsky 1995) 
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2f + 3^+1 

and the Einstein-Poisson equations for the metric or potential perturbations 

A;2$ = 4TTGa'^PT[6T + 3-(l + Wt)^], 

a k 



2^ 

{i>2) 



(A-5) 



if A; ^ ^—K. Here wt = Pt/pt where the subscript T denotes the total matter including all particle species 
and the total anisotropic stress is related to the radiation quadrupoles as 



nT = ^[02+iV2], 




(A-6) 



with N2 as the neutrino temperature quadrupole. Notice that the baryon continuity equation can be rewritten 
as Sb = —k{Vb — Oi) + 30o since dilation effects on the photon temperature and baryon density fluctuations 
are analogous. This represent adiabatic evolution if 14 = 61. 

A. 2. Acoustic Dispersion Relation 

Let us derive the dispersion relation for acoustic oscillations in the tight coupling limit. Consider first 
the efl'ects of polarization. Since only the polarization monopole Qq and quadrupole 0^ feed back into the 
temperature fluctuations, we may immediately expand the polarization hierarchy in to obtain 



e? = = ^62 



which simplifies the temperature quadrupole evolution of (A-1) to 



e2 = k 



r/202 



(A-7) 



(A-8) 



where /2 = |. Other approximations commonly used arc /2 = for impolarized radiation (Chibisov 1972) 
and /2 = 1 for further neglecting the angular dependence of Compton scattering (Weinberg 1972, Peebles 
1980, HSa). We keep the factor /2 imphcit so that the separation of efl'ects can be read directly off the final 
results. Expanding the quadrupole equation (A-8) to first order in f ~^ we obtain 

-i2, 



02=f-V2"'ofc01- 



(A-9) 
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Fig. 7. Photon difFusion scale. The photon diffusion scale grows rapidly near last scattering 
due to the increasing mean free path of the photons but remains well under the horizon scale 
= {a/a)\a, at last scattering. The small difference between a* and aa is sufficient to 
cause a significant difference in the effective damping if ili,h^ differs substantially from 
the crossover point 0.03. The inclusion of the angular dependence of Compton scattering 
enhances damping by a small factor /2 = 9/10 as does the further inclusion of polarization 
/2 = 3/4. 

A second-order expansion for the quadrupole is not necessary since its effect on the fluid equations through 
6i is already of first order in f"^. 

On the other hand, a second order expansion of the baryon Euler equation is necessary. Let us try a 
solution of the form Qi oc expz J ujd/q and ignore variations on the expansion time scale d/a in comparison 
with those at the oscillation frequency uj. The electron velocity, obtained by iteration of the Euler equation, 
is to second order 

Vb = Oi - f-^R[iujei - /c*] - T-^R^uj^Qi. (A-10) 
Substituting this into the photon dipole equation (A-1) and eliminating the zeroth order term yields 

4 



iuj{l + i?)ei = k[Qo + (1 + i?)*] - f-^iiVei - —f-^f-^k'^Q^. 
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(A-11) 



This suggests that we try a solution of the form ©o + (1 + i?)^' oc expi J udrj. Employing the monopolc 
equation of (A-1) and again assuming that variations at the oscillation frequency are sufficiently rapid that 
changes in i?, and \1/ can be neglected, we obtain 



(A-12) 



With the first order relation u) = k + R), the solution to the resultant quadratic equation is 



(l + i?)2^5-'2 l + R 



(A-13) 



18 



A n^=o.2; 



0.5 



-0.5 




10- 



10- 



10- 




(b) Spectrum z = 2000 



10-2 



10-1 
k (Mpc-') 



Fig. 8. Baryon drag effect in adiabatic CDM models, (a) Baryons cause a drag effect on 

the photons leading to a temperature enhancement of —R^ = l-R^*! inside potential wells 
which shifts the zero point of the oscillation (short dashed lines), (b) This contribution 
yields alternating peak heights in the rms and is also retained after diffusion damping. Here 
numerical results are displayed. 

Thus to second order acoustic oscillations are damped as exp[— (fc/fc£))^] with the damping length (Weinberg 
1972, Kaiser 1983, Bond 1995) 



li?2+4/-i(i + ij)/5 



(1 + Rf 



(A-14) 



If i? < 1, polarization and the angular dependence of Compton scattering enhances damping through the 
generation of viscosity (see Fig. 7). Viscosity is related to photon diffusion because the quadrupole and 
higher moments are generated as photons from regions of different temperatures meet. 

A. 3. Baryon Drag and the Adiabatic Invariant 

The temperature perturbation oscillates around 6o -I- (1 -I- i?)vE' = due to the baryon drag effect. This 
can be more easily understood by examining the first order equation from which it originates, 



(i + i?)eo + yeo 



-{l + R)^, 



(A-15) 



ignoring slow changes in i?, $ and ^ from the expansion. Notice that mcff = {\ -\- R) plays the role of the 
effective mass of the oscillator and the gravitational potential provides the effective acceleration through 
infall. The photon pressure acts as the restoring force and is independent of the baryon content R. 
Equation (A-15) has the immediate solution 



1 . 



Oo = [60(0, fc) + (1 + i?)*] cos(a;?7) + -60(0, k) sm{ojr{) - (1 + R)^ 



(A-16) 



where the frequency uj = k/\/Wn^ = k/ •\/3(l -|- R) in agreement with the first order dispersion relation of 
(A-13). This solution describes an oscillator whose zero point has been displaced by — mes^' = —{1+R)'^ due 
to the gravitational force. The photons, although massless, suffer infall effects due to gravitational blueshift. 
Since this is exactly cancelled as the photons stream out of the potential wells, we can consider Go + ^' as 
the effective temperature perturbation. Thus this part of the zero point shift has no net effect. However, the 
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baryons also contribute to the effective mass of the fluid. Since the photons and baryons are tightly coupled, 
baryonic infall drags the photons into potential wells and contributes —R^ to the displacement. Notice 
that baryon drag also increases the amplitude of the cosine oscillation since the initial conditions 8(0, fc) 
represent a greater displacement from the zero point. Thus baryon drag accounts for both the alternating 
peak heights of the acoustic oscillations and their enhancement with Q^h^ (IfSa). After diffusion damping 
has eliminated the oscillations themselves, the zero point shift — iJ^* remains. Of course, for adiabatic BDM 
models ^ is also reduced to zero by the diffusion. 

In reality, the variation of oscillator parameters, such as the effective mass and gravitational force, on 
the expansion time scale cannot be ignored over many periods of the oscillation. Prom equations (A-1) and 
(A-3), the full first order equation is 

^(1 + R)eo + y eo = -y (1 + i?)M/ - ^(1 + (A-17) 

where the addition of the space curvature term comes from the dilation effect Go = — ^. Notice that the left 
hand side is precisely the equation of an oscillator with a time- varying mass rrieff = (l+i?). The homogeneous 
equation can be solved by employing the fact that variations over a single period of the oscillation are small. 
The adiabatic invariant associated with an oscillator is given by the energy E = ^rrieai^^-A? over the frequency 

to. The amplitude therefore scales as A oc cj-'/^ oc (1 + This yields fundamental solutions of the form 

(1 + i?)-i/'^exp(±i / iodri) (Peebles & Yu 1970). The phase integral can be performed analytically 

j udri — ^ J ^sdr] = kvs, 

where rg is the sound horizon given in equation (1). 
A. 4. Summary 

Several points are worth emphasising: 

(1) The flrst order dispersion relation for acoustic oscillations is a; = kcs = k/^3{l + R). 

(2) Slow changes in the baryonic contribution to the effective mass cause the temperature oscillation to decay 

as (l + i?)-i/4. 

(3) The oscillation phase is related to the sound horizon Vs = J Cgdrj by f udrj = krg. 

(4) Photon diffusion alters the dispersion relation and leads to exponential damping. 

(5) The damping length increases roughly as Ad ^ k~j^ ~ oi' the geometric mean of the conformal 

time and the Compton mean free path as one expects of a random walk. 

(6) If it is well under the horizon, the damping length surpasses the wavelength when f/k = kr]^loT when 

the optical depth through a wavelength is still high. 

(7) The angular dependence of Compton scattering and polarization increases the damping length when 

R<1. 

(8) The zero point of oscillations in 6o + * is — ii* due to baryon drag and remains as a temperature shift 

after diffusion damping. 
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Appendix B. Horizon Crossing 

B.l. Photon-Baryon System 

The amplitude of the acoustic oscillation is determined by the growth of fluctuations before and in par- 
ticular during the epoch when the scale crosses the horizon. Since the second-order tight-coupling expansion 
just yields a multiplicative diffusion damping factor, let us obtain the first order solution which we denote 
by overhats. The full solution can be constructed through the relations 

So + * = [Oo + (1 + i?)*] exp[-(/e/fcz3)2] - m, 
61 =eiexp[-(fc/fcD)2], 

(B-1) 

56 = 3eo + 5(0,fc), 
kVb = fcSi = -Go - ^, 

where recall S{0,k) is the initial entropy perturbation in the photon-baryon system 5(0, fc) = Sb{0,k) — 
|(5-y(0, fc). Near or above the horizon at ry*, where the sources from potential growth and decay drive the 
oscillator, the full formalism of HSa is needed to follow the fluctuations accurately. However, we are only 
interested in small scale fluctuations which enter the horizon well before last scattering. Since radiation 
pressure prevents the growth of fluctuations in the radiation dominated epoch, the gravitational potentials 
decay away after horizon crossing. On small scales, the acoustic oscillations therefore experience a boost at 
horizon crossing and thereafter settle into pure modes Oo('7, k) = (l+i?)^^/^[CA(fc) cos(krs)+Ci{k) sin(fcrs)]. 
Equation (A-17) yields exact solutions for C'a and C/ for the fundamental adiabatic and isocurvaturc modes 
of the fluctuation if anisotropic stress IIt is ignored so that \1/ = — The adiabatic mode arises from an 
initial curvature perturbation $(0, k), whereas the isocurvature mode from an initial entropy perturbation 
S{0,k). The solutions are (Kodama & Sasaki 1986, HSb) 

3 

lim CA(fc) = -$(0, fc), adiabatic 

lim C/ffc) = — -^S{0,k), isocurvature 

n^o 4 k 

where keq = (20o-ffo/'*eq)^^^ = 9.67 x 10~^Cloh'^{l - /l/)^/^6^7 Mpc~^ is the wavenumber that passes the 
horizon at equality. Recall that fu = pi,/{p^+pi,). The two modes stimulate pure cosine and sine modes since 
the gravitational forcing function yields near resonant driving with the phase fixed by $(0, fc) = constant and 
$(0, fc) = respectively (HSb, Hu & White 1996). The amplitude of the fiuctuations is easy to imderstand 
qualitatively. For the adiabatic case, 0o(O, fc) = i<I'(0, fc). If photon streaming is ignored G) = — the 
dilation eflfect would raise the amplitude to ©o(?y, fc) = \^{^, k) - $(??, fc) -|- $(0, fc) = |3>(0, fc). Note that a 
decaying potential boosts the acoustic amplitude due to the gravitational forcing effect. A similar analysis 
for the isocurvature mode accounting for potential growth outside the horizon explains the isocurvature 
amplitude (HSb). 

The effect of anisotropic stress can be considered as a perturbation (HSa). The dominant term comes 
from the neutrino quadrupole since photon anisotropies are damped exponentially with optical depth before 
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last scattering. The order of magnitude can be simply read off the initial conditions for the growing mode 
of the perturbation, 

1 1 / 2 

60(0, k) = --*(0, k) = - [1 + -f^) $(0, k), adiabatic (B-3) 
or 60(0, k) = *(0, k) = $(0, k) = and 



0o(O, k) = ^fceg (1 + ] S{0, k), 



$(0, k) = ^keg (^1 + S{0, k), isocurvature (B-4) 



^{0,k) = -^keq (l--^f.]S{Q,k). 



16 ^ V 15' 

Indeed, we find that the adiabatic amplitude is well approximated by 



Ca(A:) = - 1^1 + -/^j $(0,A;), adiabatic (B-5) 

and likewise 

Ci{k) = - ( 1 - ^fA S{0, k), isocurvature (B-6) 

4 fc \ 15 / 

for the isocurvature mode. Explicitly, the first order acoustic solution is 

60 + * ~ 60 = (1 + R)-^'^[Ca cos(fcr«) + Ci sin(fcr«)], 

61 = H = -\/3(l + R)-^/^[Ca sm{krs) + Ci cos(fcrs)]. 



(B-7) 



In Fig. 1, we display an adiabatic example. In this case, it is also useful to compare the acoustic amplitude 
to the Sachs-Wolfe effect (Sachs & Wolfe 1967) \Q + ^rms{m,k) = \^{rio,k). With the relation (see HSa 
eq. A18), 



2 



*('?o,^) = -Y^(i + Y^/.) (1 + 5M ^(O'fc). (B-8) 



valid at large scales A; -C k^q, the relative amplitude becomes 
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|3CA(fc)/*(?7o,A;)| =5(^1 + -/,J , (B-9) 

and represents a significant boost. 
B.2. CDM Component 

The evolution of the cold dark matter fiuctuations in the presence of acoustic oscillations is also interest- 
ing and relevant for determining the small scale behavior of the matter transfer function (sec Appendix D) . 
The cold dark matter evolution equations arc of the same form as the baryon continuity and Euler equations 
save for the absence of coupling to the photons, 

5^ + -Sc = -kH - 3$. (B-10) 
a 
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In the radiation dominated epoch, the metric terms on the right hand side are dominated by perturbations 
in the radiation and may be considered as external driving forces. The homogeneous equation has two 
fundamental solutions 5c oc {In a, 1}. The particular solution is constructed via Green's method 

(5c-Cilna + C2+ P [In a' - In a]^{k^^ + '■i<i>)dri' . (B-11) 
Jo a 

Adiabatic initial conditions require Ci = and C2 = 30o(0, fc). Thus 6c remains constant outside the 
horizon and then gets a kick from infall and dilation that generates a logarithmic growing mode. Since the 
behavior of the potentials is self-similar in k, i.e. they are constant outside the horizon and decay to zero as 
inside of it, their effect on Sc is the same for all k. Once the potentials have decayed to zero, Sc settles 
into the logarithmic growing mode as 

4^/i$(0,A;)ln(/2— ) , aH<^a^aeq (B-12) 



an 



where horizon crossing occurs at 



Off ^ 1 + ^1 + 8{k/keq)^ 

CLeq ^{k/keqY (B-13) 

keg I, ^ ), 

_ — — • « » f^eq 

By numerical calculation of the integrals in (B-11), we obtain 

h = 9.11(1 + 0.128/, + 0.029/^), 

(B-14) 

72 = 0.594(1 - 0.631/, + 0.284/2), 

valid at the 1% or better level for the full range < /, < 1. As we shall see in Appendix D, this solution 
can be joined onto the growing mode in the matter dominated epoch to describe the full time evolution of 
the CDM fluctuations. 

Appendix C. Decoupling 

C.l Photon Decoupling 

The tight coupling approximation is strictly valid only well before decoupling. However, the acoustic 
modes may be simply joined onto the free-streaming solutions once diffusion damping near decoupling has 
been taken into account. The full Boltzmann hierarchy has the formal solution (see HSa eq. 11) 

[e + *](r/o,fc,M)= / dr?[(eo + *-iMH)f-$ + 4']e-V=''(''-''«), (C-1) 
Jo 

where fc/i = k • 7 and curvature has been neglected. The terms in parenthesis contribute at last scattering 
due to weighting by the visibility function V-y = fe^'^ and the ISW metric terms play a role between last 
scattering and the present. This formal solution is made practical by replacing the sources Qq and Vi, with 
their acoustic solution at last scattering. 
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It may seem that employing the tight coupUng solution through decoupling would lead to erroneous 
results. In particular, the damping approximation should break down if the optical depth through a wave- 
length drops below unity. However, let us examine the tight coupling solution more carefully. Equation 
(B-1) implies 

rno \ (C-2) 

Jo 

Although baryon drag effects such as acoustic displacement and enhancement are already encorporated in 
the acoustic solution 9o [see equation (A-16)], the residual term appears beneath the diffusion scale. 
Let us ignore this for the moment. The effective visibility for the acoustic oscillations is given by 

= V^e-['=/'=°(^)l'. (C-3) 

This function is plotted for a given model in Fig. 2. Unlike V^,, V-y is exponentially damped at late times 
by the growing diffusion length and thus peaks at earlier times. Note that a 10% shift in redshift represents 
a factor of three in optical depth near last scattering. Thus the region where we expect the approximation 
to break down is given little weight in the integral. More specifically, the exponential damping insures that 
most contributions come from before the epoch at which the diffusion length surpasses the wavelength. As 
we have seen in Appendix A, the optical depth through a wavelength is high at this time and justifies the 
tight coupling expansion. 

The damping of acoustic modes through last scattering can in general occur due to two different mech- 
anisms working in equation (C-2): diffusion and cancellation. However the effects are of greatly unequal 
magnitude. Cancellation occurs since on small scales many wavelengths of the perturbation span the Comp- 
ton visibility function. Photons that last scattered at the crests of the perturbation will have their effect 
cancelled by those that scattered at the troughs. Mathematically this occurs in equation (C-2) because the 
oscillating plane wave is integrated over the visibility function. Cancellation leads to a power law damping of 
fluctuations as the scale decreases below the width of the visibility function. However, in the case of diffusion 
damped acoustic contributions, it is not the width of the Compton visibility function V-^ that is relevant but 
rather the acoustic visibility function V-y. As one goes to smaller and smaller scales (high fc), the width of 
this function decreases as well. Thus even at high k the cancellation regime is never fully reached and one 
may approximate the integral (C-2) by replacing V-^ by a delta function, i.e. 

[6 + *](%, ^ [Oo + * - inVb\{v*,k)e"'^^'"-''°'>V^{k), (C-4) 

where 

V^{k) = / dr^% = / rir/V,,e-['=^/'=°('')I . (C-5) 
Jo J a 

The observable anisotropy follows by decomposing equation (C-4) into Lcgendre moments O = ^(— i)'^0/^_P^(/i) 
and summing over fc- modes = ^ J k^\Qi{rji:),k) / {2t+\)\'^dh\k (HSa). Decoupling thus increases the effec- 
tive diffusion damping length due to the corresponding increase in the mean free path of the photons. The 
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Fig. 9. The residual baryon drag effect after last scattering in an adiabatic CDM model. On 

scales under the width of the visibility function, cancellation between contributions which 
came from potential wells and hills at last scattering damps fluctuations from the baryon 
drag effect. Note that cancellation damping is weak and scales as {kr]*)~^/'^ in contrast 
to the exponential diffusion damping. Projection relates the rms fluctuation in (a) to the 
anisotropy power spectrum in (b). 

result is a near exponential damping with scale that completely overwhelms the small residual cancellation 
damping. 

Cancellation damping docs occur for the baryon drag effect —R^ in equation (C -2) which remains after 
diffusion damping. The amplitude of the resultant fluctuations can be estimated by noting that 

rvo 

[6 + *](?7o, k,n) = - / dr]V^me''"'^''-''°^ (C-6) 
Jo 

is approximately a Fourier transform. Employing Parceval's theorem, we obtain 

1 

~ J r |i?*V^|2d77 (C-7) 
1^ Jo 

Thus the contribution to the rms is suppressed by roughly due to cancellation. Note that the resid- 

ual baryon drag effect only appears in models where ^ itself is not damped away by diffusion, e.g. adiabatic 
CDM and isocurvature BDM models. In Fig. 9a, we display this effect. The amplitude of the effect is slightly 
overestimated since baryon drag weakens through last scattering. Its effect on the anisotropy is shown in 
Fig. 9b and can be obtained analytically in a computationally simple matter through the approximations of 
Hu & White (1995). Since it is unlikely to be observable due to effects in the foreground of last scattering, 
we omit a detailed calculation here. 
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C.2 Baryon Decoupling 

The formal solution to the baryon Euler equation (A-3) is 

aVb^ f dri' a[fdei+k^]e-^\ (C 8) 

Jo 

where recall fd = f/R and quantities in the integrand are evaluated at rj' . The two source terms are the 
Compton drag effect and infall into potential wells. The former forces the baryon velocity to follow the 
photon dipolc (velocity) at high drag optical depth r^. The presence of the scale factor a in the equation 
represents the fact that baryon velocities decay as a^"^ in the absence of sources. Since fde^^'' is very nearly 
a delta function with respect to variations on the expansion time scale, this equation is conceptually identical 
to its photon analogue (C-1) with the replacement r — > fd. The plane wave factor is absent for the baryons 
since their particle velocities are low and the streaming can be neglected in comparison to the wavelength. 
Drag depth unity Td{zd) = 1 marks the transition between the drag and infall epochs. For > 0.03, the 
drag epoch precedes last scattering Zd > z* assuming standard recombination. If the universe is reionized 
after recombination to some ionization level .Xg, Zd = 263{Cloh'^)^'"^Xe^^^ {1 — Yp/2)^^^^'02 9,nd the drag 
epoch significantly precedes last scattering in most scenarios. t Here Yp is the helium mass fraction Yp ~ 0.23 
By analogy to the photon case, it is useful to define the drag visibility function 

suitably normalized to have unity area. Diffusion damping modifies the acoustic visibility function as 

V(, = Vfee-l'=/'=° (C-10) 

Thus, we expect that immediately after the drag epoch the baryon velocity and density perturbations are 
approximately 

Vb{m,k) = ei{vd,k)Vb{k), 

(C-11) 

h{Vd,k) = db{vd,k)'Db{k), 

where 

Vbik) = / drjVb = / dryV6e-['=/'=^('')J , (C-12) 
Jo Jo 

and recall that ©i and Sb were given in (B-7). 

Appendix D. Matter Evolution 

After horizon crossing but before the end of the drag epoch, baryons follow the acoustic solution of 
(C-11) and the CDM follows their own pressureless evolution. After the drag epoch, the baryon evolution 
equation (A-3) is identical to the cold dark matter and their joint evolution can be expressed in terms of 
fluctuations in the total non-relativistic matter density dm- Thus, the solution for the time evolution of the 
matter fluctuations requires knowledge of both the baryon and CDM perturbations at the drag epoch. The 
bar>-()ni(: contribution was olitaincd in Appendix C. Let us now evaluate the CDM contributions. 



^ This differs from the treatment of (HSb) where Zd was deflned to be the epoch when the perturbation 
joined the growing mode of pressureless linear theory. The presence of a decaying mode lowers this redshift 
by a factor of | [see Appendix D, eqn. (D-19)]. 
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Fig. 10. CDM evolution in the Compton drag epoch. If baryons contribute a significant 
fraction of the total matter density, CDM growth will be slowed between equality and the 
drag epoch. Held by Compton drag, the baryons do not contribute their self-gravity. For 
the numerical results, we choose a model that never recombined so that ^ a^q. 

D.l Exact CDM Solutions 

The evolution of CDM fluctuations is described by equation (B-10). Since the curvature or A terms are 
negligible before the drag epoch, this equation can be rewritten in terms of the equality-normalized scale 
factor y = a/ueq as 



^ , (2 + 3;/) d ^ _ 



(D-1) 



2y{l + y) Qo 

Here we have assumed that the radiation contributions to the gravitational potential have decayed to zero 
well after horizon crossing. In typical adiabatic models, the CDM contribution usually dominates the non- 

relativistic matter. Consider first the limit of negligible baryon fraction, fib/Oo <C 1. In this case, the 
matching condition at the drag epoch becomes trivial since the baryons have no effect on the CDM. If 
Qc = ^^Oi equation (D-1) has same solution before and after the drag epoch (see Peebles 1980 eqns. 12.5, 
12.9), 

_ 2 

■y, 



U2 



1 K 

= ^(2 + 3y)ln 



f(l+y)^/^ 



_(l+y)l/2_l_ 

before curvature or A domination. Matching to the radiation dominated solution (B-12), we obtain 



(D-3) 



for k, 



eq- 



Let us now solve the equation (D-1) for the case that the contribution of baryon is not negligible. The 
two independent solutions are given in exact form through Gauss' hypergeometric function F by 



Ui = il+ y^'Fiai, ai + ^,2ai + ^) , 



(D-4) 
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where i = 1, 2 and 



Cti 



1 ± VI + 24nc/flo 



(D-5) 



with — and + for i = 1 and 2, respectively. Note that liniy^ooUi = y~°'* . Thus the main effect of the 
baryons is to slow the power law growth of CDM after equality. 

It is easy show that these solutions are identical to equations (D-2) for fic = Oq. They also take on 
elementary forms for two other special cases: Q,c = 0, 



Ut =1, 



C/2 =-ln 



(l+t/)V2_l 



(D-6) 



and Oc = k^Oi 



(l+?;)i/2 + l 



(D-7) 

.(l+?y)i/2_l_' 

In order to map the solution for the radiation dominated limit (B 12) onto amplitudes of Ui and U2, we 
have to take the limit as y ^ of (D-4). By using a linear transformation of the hypergeometric function 
(see e.g. [1] eq. 15.3.9), we find 



lim Ui 



T{2a, + 1/2) 
r(a,)r(a, + 1/2) 



- Iny + 2V(1) - V(ai) - V'l^i + 1/2)] 



(D-8) 



where r(x) and ijj{x) = r'(x)/r(a;) are gamma and digamma functions respectively. Matching to the 
radiation dominated solution (B-12) yields 



Sciv, k) = 7i$(0, k) [AiC/i(?7) + ^21/2(77)] 



(D-9) 



where 



Ai=- 



r(ai)r(ai + 1/2) 



r(2ai + 1/2) [^{ai) + V(ai + 1/2) - ^(a2) - ^(^2 + 1/2)] 



In ( ) + 2V(1) - V(a2) - V'(a2 + 1/2) 



(D-10) 



A2 is obtained by replacing the subscripts 1 ^ 2 of Since fib/rio ^ 1) it is useful to approximate the 
coefficients with a series expansion, 

B2, (D-11) 



Ai=Bi In 



aH 



with 



Bi = -[l- O.568(06/f^o) + 0.094(^6/^0)^ + 0[{Qb/^of]] , 

B2 = ^(ln4 - 3) [1 - 1.156(17b/f}o) + 0.149(f}fe/f7o)' - 0.074(nb/17o)3 + O[(0b/l^o)^]] , 



valid at the percent level for Ob/f^o < 1/2, and 



r(a2)r[a2 + 1/2] 
r[2a2 + 1/2] 



r(2ai + 1/2) 
r(ai)r(ai + 1/2) 
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Fig. 11. Velocity overshoot effect. Below the horizon at the drag epoch krjd ^ 1, the 
acoustic velocity at Zd dominates the growing mode and hence the final transfer function. 
Neax the horizon, the acoustic density becomes comparable and shifts the zero points of the 
oscillation. 

For the three special cases, we can describe Sc by elementary functions. If Oc = f^Oi is given by equation 
(D-2), whereas 

Oc = , 



Sc{v,k) =/i$(0, k) 
6c{r],k) =/i$(0,/c) 



In (4/2^ ]Ui-2U2 

an 



ln{Ahe-'^)U,-'-U, 
aH J 3 



1 



(D-14) 



• 



In Fig. 10, we show the time evolution of 5c before the drag epoch for several different values of Q.b/0^c- 
Numerical results in this figure arc for fully reionized models so that the drag epoch ends well after equality 
unlike other examples in this paper. 

D.2 Matter Transfer Function 

With the baryon and CDM fluctuations at the drag epoch from equations (C 11) and (D 9) respectively, 
we can now solve for the evolution to the present. After the drag epoch, baryons behave dynamically as 
CDM and the combined non-relativistic matter fluctuations 



(D-15) 



follow the growing and decaying solutions for 5m (Peebles 1980) 

2 



Do 



15 



(2 + 3y)ln 



_(l+y)l/2_l_ 



45 



+ 



1/2 



(D-16) 
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before curvature or A domination. To account for effects from curvature and A at a » a^q ov y ^ 1, one 
simply needs to replace D, where 



D{a) = ^Oo5(o) / /° ..3 , 



(D-17) 



a ^Oo + a ^(1 - - ^^a) + f^A, 
is the growing mode of radiationlcss linear theory normalized to equal a at early times. 
By matching the fluctuations at the drag epoch, we obtain for the growing mode 

D2 ^ , ^ D2 



(D-18) 



D1D2 - D1D2 D1D2 - D-i_D2 

and similarly for the decaying mode. The combined time evolution is plotted in Fig. 4 for a BDM model 
with 5m = 5i). If Z4 <^ Zeq, equation (D-18) reduces to the familiar form 



5m{r],k) 



a 
a-d 



3 1 
<5m(%, k) - -{kr]d)Vm{Vd, k) 



(D-19) 



_5 ^ 5' 

Notice that on scales much less than the horizon at the drag epoch, the velocity at % dominates the growing 
mode if the two values are comparable at r]d (see Fig. 11). This "velocity overshoot" effect occurs since the 
peculiar velocity moves the matter and creates density fluctuations kinematically. Expansion drag on the 
velocity eliminates it in an expansion time 77^ and thus causality prevents this effect from generating density 
fluctuation above the horizon at the drag epoch krjd <C 1. 

It is conventional to recast the evolutionary effects in terms of a transfer function. As with equation 
(D-15), we can break the present day transfer function up into a baryonic and cold dark matter contribution 
at the drag epoch 

o,. / o,_ \ 

(D-20) 



"0 



1 - ) Uk). 



It should be kept in mind that Ti, and Tc do not represent the respective transfer functions today. Let us 
consider the adiabatic transfer function. Here one expresses the evolution of small scale fluctuations in terms 
of those at large scales, i.e. \5T{'r}o,k)\'^ oc T'^{k)P{k) with normalization limfe^o ^'(A;) = 1 and initial power 
spectrum P{k) oc |fc^$(0, k)\^. The large scale solution is given by (see HSa eqn. A-15) 



lim 6t{t], k) 
fe— >o 



|(-^) A(7?)$(0,fc). 



5 Vfc, 



The acoustic contribution of the baryons is therefore 

2 



("^') '^b{k){l+R)-^''^ (^oskrs - ^kcssinkr^ Gi 



(D-21) 



(D-22) 



This equation is compared with numerical results for the adiabatic transfer function in Fig. 5a. 

For isocurvature BDM models, T = and it is conventional to define it such that |^t(?7o,^)|^ oc 
T^(/c)|5(0, fc)|^ with normalization limfc_>oo T'(fc) = 1- From equation (B-6), the small scale tail of the 
transfer function becomes 



T{k) = 1 



3V6 



1 



15 



-^Vt,{k){l + R)-^'^ 



sin kvs 



£1 



-kcs cos krt 



(D-23) 
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This function is plotted in Fig. 5b. 

In the 0(,/0o ~* limit, the CDM contributions can be expressed in terms of elementary functions as 



Inl.Sg 



(D-24) 



where q = {k/Mpc ^)Q2.7i^ah^) ^- This should be compared with the high k tail of the standard fitting 
function to the numerical results (BBKS), 

rBBKs(«) = ^''^^34^'^^ ^^ ^ ^■^^'^ ^ ^ ^ (6.71g)4]-i/4, (d-25) 

i.e. limg^oo T'C?) = ln(2.34g)/15.7g'^ which differs by 10% from the analytic prediction at small scales. In 
fact, since the fitting formula was designed to fit intermediate scales, equation (D-24) is more accurate at 
extremely small scales. 

If the baryon fraction is non- negligible, the contribution is expressed in terms of hypergeometric functions 
through equation (D-4), 

Te = /i (^1 + ^/.) (^1 + I (^^) {Gi[AiUi + A2U2] - G2[Ai{7i + A2U2]} |^^^^. (D-26) 

Though exact, this expression is rather complicated. It is useful and instructive to seek a simple scaling 
relation for this form. Notice that for all cases the scale dependence of the transfer function may be written 
as 



Tc^a-^, (D-27) 



lnl.8/3g 
14.2g2 

for fi, = 0.405. Here a and /? are functions of Oq^^ and flh/^o- Note that as g — > 00 the modification due 

to (3 becomes insignificant. A very accurate fit to both a and /3 is given in Appendix E eqns. (E-12) (E-13). 
The numerical, analytic and fitted analytic results are compared with the empirical scalings of (Peacock & 
Dodds 1994, Sugiyama 1995) in Fig. 12. The analytic calculation is essentially exact while the fitted analytic 
form works to 1% accuracy. Notice that in this extreme case we have significantly improved upon previous 
results. 

Equation (D-27) breaks down for intermediate to large scales. The CDM contribution can be approxi- 
mately scaled from the BBKS form as 

Tc{k) ~ TBBKsiq), 

This expression is employed in Fig. 5 with the coefficient 6.71 in (D-25) replaced by 6.71(14.2/15.7) = 6.07 
to match the analytic small scale tail. Notice that at the largest scales, this expression underestimates the 
matter transfer function. This is because baryon contributions must be properly included. Although the 
limiting form limfe_>o = 1 is simple, the behavior near the horizon scale at Zd is not. Since this region is 
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Fig. 12. Adiabatic CDM transfer function in a high flb/^lo = 2/3 case. The analytic 
solution is essentially exact in the small scale limit. Simple fits based on the BBKS form 
can cause large errors at the small scale: PD (Peacock & Dodds 1994) and S (Sugiyama 
1995). The fitting function developed here [see equation (D-27)] works at the 1% level even 
for this extreme case. 

not the main focus of this work, we do not attempt to describe this analytically. If the baryonic oscillations 
are small or smoothed over, an approximate patch is given by 



T{k) ~ TBBKs(g) 



-1/2 



1 - ^ ) a-V2(Oo/.2)-iei,„ 



(D-29) 



which extends the Peacock & Dodds (1994) approach to high Qt/^o- 

Finally, we can express the ratio of the acoustic peak heights to the CDM tail with equations (D-20), 
(D-22) and (D-24), 



V3 



^A:G2(%)[1 +i?(%)]-3/4l?,(A;) h + 



In ( 4/26-3^ 



lima.^o Tc 



(D-30) 



if the velocity overshoot effect dominates the acoustic contributions. We can simplify this expression by 
noting that G2(??d) — |(0oilo/'^eg)~^^^(l + ^d/aeg)"^/^. Furthermore, the function kVb peaks at roughly 
O.Sfcs with an amplitude of OAks- With the scaling of equation (D-27), the peak relative amplitude of the 
acoustic oscillation is approximately 
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(D-31) 



The logarithmic term is roughly unity and may be dropped for estimation purposes [c/. equation (9)]. Equa- 
tion (D-31) roughly quantifies the prominence of the acoustic oscillations in a CDM model. For best accuracy 
however, the solutions (D-22) and (D-26) for the baryons and CDM respectively should be employed. 
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Appendix E: Recombination Fitting Formulae 

Rather than recalculate the atomic physics of recombination each time one needs to consider effects at 
the last scattering and drag epochs, it is convenient to have accurate fitting formulae that encorporate the 
ionization history. In general, all quantities associated with the ionization history must bo functions of iloh^ 
and rifcft,^ alone once the CMB temperature To = 2.726K (Mather et al. 1994), neutrino fraction fi, = 0.405, 
and helium fraction Yp ~ 0.23 are fixed. Fitting functions in this Appendix are designed to be valid at the 
percent level for an extended range of parameter space, 0.0025 < flbh'^ ^ 0.25 and 0.025 < floh"^ < 0.64 and 
consequently appear rather complicated. We employ a recombination calculation based on the improvements 
discussed in Hu et al. (1995). 

The last scattering epoch is a very weak function of parameters and is given by 

= 1048[1 + 0.00124(O6/i2)-0-738][i ^ g-^^^oh'^y^ 
gi = 0.0783(06/1^)"°-^^^ [1 + 39.5{Qbhy-^^^]-\ (E-1) 
92 = 0.560[1 + 21.1{Qi,h'^y-^^]-\ 
The drag epoch ends at a related redshift which depends somewhat more strongly on the parameters 

h = o.3i3(fio/i')-°-^'^[i + om7{noh''f-^'% (^"2) 

62 = 0.238(r2o/i^)°-^^^. 

The two are approximately equal if Qbh"^ ~ 0.03. 

The diffusion damping envelope can be approximated through the first decade of damping by the form 

P^(A:) ~ e-['=/'=°^l'"\ (E-3) 

where the effective damping scale kn-y has simple asymptotic scaling, 

kpj ^ / Fi, Qbh^ > 0.1 , , 

Mpc-^ I {i^bh^r'F2, nbh? « 0.1 ^> 

with 

Fi = 0.293(Oo/i^)"-^'^^[l + (25.10o/i2)"'"'^% 

F2 = G.524(Oo/i2)0-505[i + (10.500/1')"°-^'^^], (E-5) 
Pi = 0.29. 

The basic scaling in the low Vlbh? limit can be understood by the Saha approximation in which the ioniza- 
tion fraction approximately scales as Xe oc (Ofc/i^)^^/^. Thus the diffusion length ^ ^ ^f^*!^ ^ 
^J'^i^b^)^^^'^ ■ Since 77* oc (Oo/i^)~^/^ in the matter dominated high Oq/i^ limit, this is approximately of 
the same form as equation (E-5). For high Ofeft,^, the corrections from an accurate treatment of the atomic 
levels becomes more important due to the high Lyman-a opacity. These two simple limits can be accurately 
joined by a rather artificial looking but highly accurate form 
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(F2/Fi)P^/p^ {Q.bh^y^ 



Pl/P2 



Mpc-1 \7r L2^ " V o-y jj i. ^5,-6) 

P2 = 2.38(fio/i^)°-'^^ 
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Prom Fig. 2, we see that this overestimates the true damping as Clbh^ — > due to a breakdown of tight 
coupHng. It is interesting to note that the full numerical results suggest that the Saha estimation of ~ 0.25 

in equations (E-5) and (E-6) is a somewhat better phenomenological fit. 

The steepness index m~^, of the diffusion damping envelope is a very weak function of cosmological 
parameters. In the limit that last scattering occurred instantaneously m-y 2. The finite width of the 
visibility function modifies this as 

m-y = 1.46(Oo/i^)° °^°^ (1 + 0.128arctan{ln[(32.8Q6/i^)"°-*^^^]}) , (E-7) 

which only varies by ~ 10% across the full range of parameter space. 
Silk damping for the baryons can likewise be approximated by 

p6(fc)~e-('=/'=^)'"', (E-8) 

with 

and the steepness index by 

(O,/,^)-0.0297(O„fe2)0.0282 

- ^-^^ 1 + (781fi./.2)-o.926 • (E-10) 

As is the case with the photons, the latter accounts for the width of the visibility function and is almost 
independent of cosmological parameters. 

Finally, by employing equation (E-2) for the drag epoch, the cumbersome analytic result for the CDM 
dra^f-contribution to the small scale transfer function T = (f26/f2o)^6 + (^c/^o)^c from (D-26) can be fit as 

ln(L8^ 



with 



-Oi,/Oo„-(Of,/Qo)-' 
''1 0^2 



a = Qi 

01 = (46.9r!o/^')"-'™[l + (32.If7o/i2)-0-532]^ (e_i2) 

02 = (12.0fioft2)0.424j^ ^ (45.0i7o/i')-°-'''], 



as the suppression factor and 



/3-l=l + 6l[(17eM)'=-l], 

bi = G.944[I + (458f^o/i^)"°'™1"\ (E-13) 
62 = (0.395^7o/^2)-°■°2*'^ 

as the correction to the logarithm. 
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Symbol 


Definition 


Equation 


e 


CMB temperature pert. 


(A-1) 


©0 


CMB monopole pert. 


(A-1) 




y^-CMB dipole pert. 


(A-1) 




Anisotropic stress pert. 


(A-5) 




Curvature perturbation 


(A-5) 




Newtonian potential 


(A 5) 




Critical fraction in i 


(1) 




Baryon T-suppression 


(E-11) 




U1U2 power law 


(D-5) 





Baryon T-log correct. 


(E 11) 


n (rit) 


Conformal time (at t) 


(A-1) 


k 


Density pert, in i 


(8) 


T 


Compton optical depth 


(A-1) 


Td 


Drag optical depth 


(C-8) 


LU 


Acoustic frequency 


(A-13) 




Photon damping factor 


(C-5) 


Vb 


Baryon damping factor 


(C 12) 




Compton visibility 


(C-1) 


v., 


CMB acoustic visibility 


(C-3) 


V,, 


Drag \"isil)ilily 


[V 9) 


Vh 


Baryon acoustic visibility 


(C-11) 


Ai A2 


Boost CDM matching cond. 


(D-10) 


Ca 


Adiabatic acoustic amp. 


(B 5) 


Cj 


Isocurvature acoustic amp. 


(B-6) 


D 


Radiationless growth 


J- 1 ) 


-Di Do 


Matter-radiation growth 


fD-16) 




Drag matching cond. 


fD-18) 


h h 


CDM boost integrals 


(B 14) 


Ui U2 


Drag CDM growth 


(D-4) 


RiRt) 


6/7 momentum (at t) 


(A-2) 


T 


Transfer function 


(D-20) 


n 


Baryon T drag contrib. 


(D-22) (D-23) 


Tc 


CDM T drag contrib. 


(D-28) 


Vi 


Velocity in i 


(A-1) 



Table 1: Commonly used symbols.^ 

t Fluid elements are i = 6, c, 7, i^, to, r for the baryons, CDM, photons, neutrinos, total non-relativistic 
matter, and total relativistic matter respectively. Special epochs include t = *,d,eq,H for last scattering, 
the drag epoch, matter-radiation equality, and horizon crossing respectively. Overdots represent conformal 
time derivatives. 
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Table 1 (cont.) 



a (at) 


Scale factor (at t) 


(A-1) 


an 


Horizon crossing a 


(B 13) 




Photon-baryon sound speed 


(A-2) 


u 


Neutrino fraction 


(B-2) 


k 


wavenumber (Lapl. eigenv.) 


(A-1) 




Diffusion wavcn. 


(A-14) 


ko-f 


CMB damping waven. 


(A-14) (E-6) 


ks 


Silk damping waven. 


(6) (E-9) 


keq 


Equality horizon waven. 


(B-2) 


k-yj 


jth photon peak waven. 


(2) 


kt)j 


jth baryon peak waven. 


(5) 


m-y 


CMB damping steepness 


(4) (E-- 7) 


ms 


Baryon damping steepness 


(6) (E-10) 


q 


Scaled wavenumber 


(D-25) 




Sound horizon 


(A-2) 


y 


a/Ueq 


(D-1) 




Last scatt. redshift 


(E-1) 


Zd 


Drag redshift 


(E -2) 


Zeq 


Equality redshift 


(1) 
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